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ABSTRACT 

As  a  port  of  the  Laboratory's  fundamental  shielding  studies  for  personnel 
shelters,  fast  neutron  dose  rotes  are  calculated  in  the  second  leg  of  an  air  duct 
through  concrete  for  neutron  energies  of  14  Mev  and  2.5  Mev.  The  colculational 
technique  is  based  on  the  albedo  concept.  Dose  rates  are  also  calculated  by  a 
Monte  Carlo  technique,  and  the  results  obtained  by  the  two  theoretical  methods 
are  compared  with  each  other  and  with  experimental  measuremerSts. 

Comparison  shows  very  good  agreement  among  these  three  independent 
determinations. 


Qualified  requesters  may  obtain  copies  of  this  report  from  DDC. 
The  Laborotory  invites  comment  on  this  report,  porticularly  on  the 
results  obtoined  by  those  who  have  applied  the  Information. 
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INTRODUCTION 


An  important  problem  in  shelter  shielding  is  the  streaming  of  fast  neutrons 
through  the  shelter  entranceway.  Tentative  experimental  measurements  have  been 
made  of  fast  neutron  dose  distributions  in  a  duct,  ^  but  a  theory  for  calculating  the 
dose  remains  to  be  developed. 

The  problem  of  fast  neutron  streaming  through  ducts  is  different  in  many 
respects  from  the  problem  of  deep  penetration  of  neutrons  within  a  medium.  In  the 
case  of  deep  p  •:netration,  the  interaction  of  the  neutron  with  nuclei  of  the  trans¬ 
porting  medium  is  important.  However,  for  the  duct  streaming  problem,  the  principal 
factor  is  the  reflection  of  neutrons  from  the  walls  of  the  duct.  For  deep  penetration, 
then,  the  scottering  medium  can  be  treated  as  homogeneous;  whereas  in  the  duct 
streaming  problem  the  material  of  which  the  walls  are  built  and  the  material  with 
which  the  duct  is  filled  must  be  treated  separately.  In  this  study  we  shall  be  con¬ 
cerned  with  air  ducts  through  concrete. 

Because  of  the  difference  in  nature  between  neutron  streaming  and  deep 
penetration,  the  Boltzmann  transport  equation,  which  describes  the  deep  penetration 
problem,  cannot  be  used  for  the  duct  streaming  problem.  Therefore,  special  methods 
have  been  devised  for  treatment  of  the  cose  of  a  cylindrical  duct.^  Since  the 
geometry  is  more  complicated  for  a  rectangular  duct  than  for  a  cylindrical  one,  the 
analytical  approach  has  usually  been  abandoned  in  favor  of  Monte  Carlo  techniques 
In  treatments  of  rectangular  ducts.  In  this  paper  an  analytical  approach,  using  the 
albedo  concept,  will  be  applied  to  rectangular  ducts. 

In  shielding  calculations  it  is  often  preferable  to  be  concerned  with  neutron 
dose  instead  of  neutron  flux,  because  biological  hazord  Is  more  readily  determined 
in  terms  of  dose.*  Moreover,  colculotional  techniques  are  frequently  simpler  when 
dose  Is  considered  instead  of  flux.  Note,  for  example,  that  neutron  flux  depends 
upon  neutron  energy  as  well  as  the  angular  and  spatial  variables,  while  dose,  being 
proportional  to  integrated  energy  flux,  depends  only  on  the  angular  and  spatial 
variables. 

In  order  to  formulate  the  equation  which  describes  the  variation  of  dose  within 
a  duct,  it  is  necessary  to  know  the  differential  dose  albedo  of  neutrons  striking  the 
wall  of  the  duct. 


*  I.e.,  absorbed,  or  rod,  dose.  Hereafter  ail  references  to  dose  will  be  to  absorbed 
dose,  unless  otherwise  specified. 


In  the  next  section,  a  formula  for  neutron  albedo  will  be  discussed.  Following 
that,  the  albedo  method  will  be  used  for  calculating  dose  within  a  duct.  Next,  dose 
determinations  by  a  Monte  Carlo  calculation  and  by  experimental  measurement  will 
be  discussed.  Finally,  the  results  obtained  by  these  three  independent  methods  will 
be  compared. 


SEMIEMPIRICAL  FORMULA  FOR  NEUTRON  DOSE  ALBEDO 

k 

Definition  of  Albedo 

The  dose  albedo  for  neutrons  is  defined  as  the  ratio  of  the  dose  reflected  from 
an  area  to  the  dose  which  is  Incident  on  the  reflecting  area.  Because  of  the  near 
equivalence  between  dose  and  the  energy  flux,  an  alternative  definition  Is  the  ratio 
of  the  energy  flux  reflected  from  an  area  to  the  energy  flux  which  is  incident  on 
the  reflecting  area.  The  definition  of  albedo  as  a  ratio,  as  given  above,  is  a  con¬ 
ventional  definition  for  total  dose  albedo. 

In  detailed  shielding  calculations,  it  is  frequently  necessary  to  know  differential 
angular  dose.  Therefore,  it  is  important  to  know  the  differential  dose  albedo.  This 
is  the  ratio  of  the  dose  reflected  per  unit  solid  angle  into  the  given  direction  from 
a  differential  area  dA  to  the  dose  incident  on  surface  area  dA,  as  shown  in  Figure  ]. 
Mathematically,  differential  dose  albedo  is  expressed  as 

/r  o  o  s  1  /'dDX 

o(E^,  e^,  e,  0)  =  (ja) 

O  O 


where 


E  =  energy  of  incident  neutron 
9^  =  polar  angle  of  incidence  of  neutron 
0  =  polar  angle  of  reflected  neutron 

<t>  =  azimuthal  angle  of  reflected  neutron 

D  =  reflected  dose  rate 

D  =  dose  rate  due  to  the  incident  neutron  beam 


o 

The  definition  given  above  corresponds  to  that  given  by  Chilton  and  Huddleston 
for  gamma-ray  differential  dose  albedo. 
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Figure  1.  Reflecfion  of  neutron  from  a  slob. 
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Frank  J .  Allen  et  al  have  given  the  results  of  Monte  Carlo  calculations  of 
neutrons  bcckscattered  from  a  semi-infinite  slab  of  concrete.  An  attempt  is  made 
below  to  devise  a  semiempirical  formula  for  the  differential  dose  albedo  of  neutrons 
which  will  fit  their  Monte  Carlo  results. 

Derivation  of  Formula 

In  the  derivation  of  the  desired  semiempirical  formula  the  following  assumptions 
arc-  made : 

1.  The  energy  dependence  and  the  spotial  dependence  of  the  function 
describing  the  reflected  neutron  flux  can  be  separated  Into  an  energy 
component  and  a  spatial  component,  as  follows; 


n  (E^^,  f,  O)  -  N  (E^)  77  (F,  O) 


o 


2.  The  scattering  of  neutrons  in  concrete  is  isotropic  in  the  laboratory  system. 
(Clearly,  this  assumption  is  not  volid  for  fast  neutrons  or  for  elastic 
scattering  of  neutrons  by  light  nuclei.  However,  it  may  be  essentially 
valid  for  neutron  scottering  in  concrete.) 

3.  The  angular  distribution  of  the  reflected  neutron  dose  is  dominated  by 
singly  scattered  neutrons.  The  scattered  neutrons  of  the  highest  energy 
will  be  those  which  have  been  singly  scattered.  Of  course,  the  multiple 
scattered  neutrons  also  contribute  to  the  reflected  dose.  However,  it  is 
assumed  that  the  contribution  from  the  singly  scattered  neutrons  is  much 
higher  than  the  cantributiort  of  the  multiple  scattered  neutrons.  (Admittedly, 
this  assumption  Is  a  weak  point  in  the  argument,  its  justification  consisting 
mainly  of  the  fact  that  correct  answers  are  obtained  to  the  problem  at  hand.) 
The  neutron  dose  can  be  separated  into  a  few  components,  such  as  a  singly 
scattered  part,  and  a  mu Itlpl 6  scattered  part.  However,  for  simplicity, 
only  the  singly  scattered  part  is  considered  in  this  paper, 

4.  Neutron  dose  is  approximately  proportional  to  neutron  energy  flux.  (This 
is  a  valid  assumption  for  fast  neutrons.  For  neutrons.jn  the  resonance  energy 
region,  the  assumption  is  poor.) 

Consider  pqw  a  neufron  scattering  within  a  differential  volume  element  dV 
which  has  a  unit  cross-sectional  area  and  thickness  dx  at  a  slant  depth  x  below  the - 
surface  of  a  concrete  slab,  as  shown  in  Figure  2.  The  probability  of  a  neutron 
traversing  without  the  interaction  the  distance  Inward  to  the  volume  element  dV 
can  be  expressed  Of, 

-Z,(E„)x 

e  ::  ^ 
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On  the  assumption  that  the  scattering  is  isotropic,  the  probability  of  a  neutron 
traversing  without  interaction  the  distance  from  dV  to  the  surface,  after  scatter,  is 
proportional  to 

-  r  (E.)  X  cos  6  sec  9 
t  I  o 

e 


where  I|,(E)  =  macroscopic  total  cross  section 

=  energy  of  bockscattering  neutrons 

Thus,  the  probability  for  a  neutron  to  traverse  a  distance  x  into  the  scattering  medium 
and  then  to  be  scottered  back  out  of  the  surface  is 


1  “I.(E  )x  -I.(E,)xcos9  sec  0 

I  V  I  O  11  o 


,  I  (E  )  e 
Air  so. 


where  is  the  macroscopic  scattering  cross  section. 

Integrating  over  x  from  zero  to  infinity,  in  order  to  have  all  the  singly  scattered 
components,  : 


CD 


I 


-Z^E  )x  -ZiE,)xc6s9  sec0 
t  o  t  1  0 


dx  /. 


Z^(E  )  +  Z.(E,)cos9  sec  9 
t  o  t  r  o 


On  the  assumption  that  ^|(Ep)  ^  above  equation  can  be  written  as 

1 


(E  ) 

4  7r  s  o 


(E  ) 

477  s'  o' 


Z.(E  )  +  Z.(E,)cos8  sec  6  Z.(E  )  (1  +  cos  6  sec  6) 
t  0  t'  I  o  t'  o  o 


A(g 


i  +  cos  9  sec  9 
o 
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where  A  (Eq)  is  an  energy-dependent  parameter,  equivalent  to  (1/4 rr) ris(Eo)/Z).(EQ)  ]  . 
But  A  (Eq)/!  +  cos  Qq  sec  6  is  just  the  probability  per  unit  solid  angle  per  unit  area 
normal  to  its  direction  of  motion  that  a  neutron  incident  at  6q  scatters  in  the 
direction  6.  Therefore,  the  differential  dose  albedc;  for  neutrons  can  be  written  os 


^  e  /  e) 

o  o 


A  (E^)  cos  6 

cos  0  +  cos  0 
o 


(1) 


It  is  not  necessary  to  specify  the  azimuthal  dependence  of  albedo  since  the  assumption 
was  made  that  the  scattering  was  isotropic  in  the  laboratory  system. 

Values  for  the  energy-dependent  parameter  A  (Eq)  were  obtained  for  each  ot 
the  incident  neutron. energies  (0.1,  0,25,  0.5,  1,  2,  3,  5,  and  14  Mev)  by  means  of 
a  least-squares  analysis  to  provide  the  best  fit  of  the  above  equation  to  the 
Monte  Carlo  data  of  Allen. ^ 

In  Allen's  data,  the  dose  reflection  factor  per  steradian  (DRF)  is  defined  as 


R..D.  j  =  1,  2,  ...  12 

ij  I  ' 


where  subscripts  I  and  j  refer  to  theTfH  energy  group  and  the  jth  angular  sector.  The 
terms  are  .  - 

6  =  polar  angle  of  reflection 

0  =  solid  angle  of  sector 

D  =  incident  dose  per  neutron 

R  =  reflection  factor 


(DRF).  = 


sec  e, 

_ 1- 

DO. 

J 


10 

r 

/ 

i  =  i 


The  angular  sectors  used  by  Allen  are  shown  in  Tables  I  and  II.  Table  I  shows  the 
angular  sectors,  average  reflection  angle,  and  solid  angle  of  each  sector  for  the 
case  of  normal  incidence.  The  average  reflection  angle  of  a  sector,  6,  is  the 
arithmetic  mean  of  the  polar  angles  of  the  two  end  points,  and  02,  of  the  sector. 
For  the  case  of  normal  incidence,  there  is  no  azimuthal  dependence  for  the  sectors. 
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Table  i.  Angular  Sector  Histogram  for  Normal  Incidence  (cos  0  ,  =  1) 

o 


Sector 

e 

!2 

Solid  Angle 

1 

8°18' 

0 

16°35.9' 

.0.2618 

2 

20°5' 

16° 35.9' 

23°33.4' 

0.2618 

■  3 

26°45' 

23°33.4' 

.  29°55.6' 

0.31416 

.  4 

32^35' 

29°55.6' 

35°14.8' 

0.31416 

..  5 

.37°36' 

.35°! 4.8'.  . 

■39°56.7': 

0.31416' 

■  6  ■  ■ 

42®5' 

39°56r7'' 

44° 13.2'  • 

■  0.31416 

46^12' 

44°13.2' 

48°ll.'l' 

0.31416 

8 

51 '^1 5' 

48°Vl.4' 

54°  18.9': 

0.52360 

9 

57°9' 

54° 18.9' 

60°0' 

0.52360 

10 

62°4r 

60°0' 

65°22,5' 

0.52360 

11 

67°57' 

65°22.5' 

70°31.7' 

0.52360 

12 

80° 16' 

70°31.7' 

90° 

2.0944 
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Table  II.  Angular  Sector  Histogram  for  Slant  Incidence 


Solid 


Sector 

0 

6| 

S2 

<^1 

<P2 

Angli 

1 

11°47' 

0 

23°33.4' 

0 

TT 

77/6 

2 

35°52' 

23°33.4' 

48°n.4' 

2it/3 

IT 

77/6 

3 

35°52' 

23°33.4' 

48°n.4' 

n/3 

2v/3 

77/6 

4 

35°52' 

23°33.4' 

48° 11.4' 

0 

7t/3 

77/6 

5 

59°22' 

48° 11.4' 

70°31.7' 

37t/4 

n 

77/6 

6 

59°22’ 

48°11.4' 

70°31.7' 

7t/2 

377/4 

77/6 

7 

59°22' 

48° 11.4' 

70°31.7' 

it/4 

77/2 

77/6 

8 

59^22' 

48° 1 1.4’ 

70°31.7' 

0 

77/4 

77/6 

9 

o 

o 

CO 

70°31.7' 

90° 

37t/4 

77 

77/6 

10 

80°  16' 

70°31.7’ 

90° 

77/2 

377/4 

- 

7/  /  0 

n 

80°16' 

70°31.7' 

■o 

o 

o 

tt/4 

77/2 

77/6 

12 

80° 16' 

70°31.7' 

90° 

0 

77/4 

77/  6 

Table  II  shows  fhe  angular  secl-or  histogram  for  the  case  of  slant  incidence. 

The  average  polar  angle,  0,  is  obtained  in  the  same  manner  as  for  Table  I.  The 
angles  0^  and  (j>2  are  the  end  points  of  the  azimuthal  angles  for  each  section. 

Since  the  dose  reflection  factor  computed  by  Allen  is  not  the  same  as  the  term 
a  (E^,  9  ,  6)  as  used  in  this  report,  his  values  must  be  converted  by  the  relationship 


a. 

J 


(DRF) 


cos  6. 

_ 1 


cos  6 


As  there  appears  to  be  no  significant  variation  of  a.  with  <p,  the  azimuthal  dependence 
of  albedo  is  neglected  throughout  this  report. 

The  least-squares  analysis  was  carried  out  on  the  Laboratory's  IBM-1620.  The 
values  thus  obtained  for  A  (Eq)  are  shown  in  Table  III  and  in  Figure  3. 

Until  now,  A  (Eq)  has  been  treated  as  an  empirical  energy-dependent  parameter. 
It  has  been  observed  that  the  curve  for  A(Eq)  as  a  function  of  Eq  can  be  fitted  by 
a  least-squares  regression  analysis  to  the  formula 


a 

A  (E  )  =  E  e 
■  o  o 


+  b  -i/F  +  cE 
Vo  o 


where  a  =  0.9719 
b  =  -2.895 
c  =  0.3417 


The  equation  for  differential  dose  albedo  can  now  be  written,  semiemperical ly. 


as 


E  cos  0 

a  (E  ,  9  ,  9)  =  - 2 

O  O  cos  0  -r  cos  0 


0.9719  -  2.895^  +  0.3417E 
V  o  o 


(2) 


Comparisons  between  the  one-parameter  formula  (Equation  1),  the  semiempirical 
formula  (Equation  2),  and  Allen's  Monte  Carlo  results  are  shown  in  the  Appendix. 
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Table  III.  Values  of  Energy-Dependent  Parameter  for  Ofie-Parameter 
Semlempirical  Formula  for  Differential  Neutron  Dose 
Albedo  on  Concrete 


A{E^) 

E  (Mev) 
0 

0.146 

0.1 

0.154 

0.25 

0.157 

0.5 

0.192 

1.0 

0.173 

2.0 

0.155 

3.0 

0.127 

5.0 

0.084 

14.0 

Figure  3.  Graph  of  values  obtained  for  A  (E^)  listed  in  Table  III. 
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CALCULATION  OF  DOSE  IN  A  DUCT 


On  the  basis  of  the  formula  for  differential  dose  albedo  just  derived,  it  is 
possible  to  develop  an  expression  governing  neutron  dose  distribution  within  a  duct.* 
First,  consider  the  most  general  problem  of  neutron  transport  through  a  medium. 
The  Boltzmann  type  transport  equation  is  regarded  as  the  complete  relation  for 
neutron,  or  radiation,  transport.  However,  the  transport  equation  cannot  be  easily 
solved  for  the  duct  problem  because  of  the  complexity  of  the  geometry  and  physics 
involved.  The  steady-state  neutron  transport  equation  can  be  written,  in  the  Integral 
form,'^  as 


n  (r,  O,  v)  =  J  J  J  J  J  v') 


K(P,  fi',  v'  -T,  0,  v)dv'dn'  +  (3) 

where  n  (r,  fi,  v)  =  mean  number  of  neutrons  per  unit  volume  per  unit 

solid  angle  per  unit  velocity  interval  at  position 
r,  with  scaler  velocity  v  In  the  direction  of  fl 

K(r',  O',  v'  -  r,  0,  v)  =  probability  per  unit^vojume  that  a  neutron  described 

by  the  coordinates  r'.  O',  and  y^is  scattered  such 
that  the  new  coordinates  are  r,  O,  and  v 

dV  =  differential  volume 


S  =  source  term  =  mean  number  of  unscattered  neutrons 
from  source  reaching  r,  with  scaler  velocity  v  In 
the  direction  O,  per  unit  volume  per  unit  solid 
angle  per  unit  velocity  interval 


*  Gamma-ray  dose  distribution  calculated  by  the  albedo  model,  using  a  method 
very  similar  to  the  one  developed  in  this  study,  gives  very  good  agreement  with 
experimental  results,  as  discussed  by  Chapman.^ 
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Changing  the  variable  from  scaler  velocity  v  to  energy  E,  the  angular  neutron 
energy  flux  can  be  written  as  En  (E,  r,  O).  Therefore,  neutron  dose,  D,  can  be 
expressed  as  the  angular  dependent  quantity: 


E 

max 

DfrTS)  =  y  f(E)En(E,  0  n)dE 

O 

where  E^ax  's  the  maximum  neutron  energy,  and  f(E)  is  the  energy-dependent 
response  function  for  neutron  dose.  In  the  monoenergetic  case,  the  above  expression 
reduces  to 


DiT,  Q)  =  f(E)En(r,  O) 


In  order  to  obtain  the  dose  transport  equation,  multiply  Equation  3  by  Ef(E) 
ond  change  the  variable  from  v  to  E.  This  gives 


n{E,r,n)Ef(E)  =  JJJdV 


Ef(E)K(E‘,  r'  -  E,  r,  0) 


n(E',’;^,  n')E'f(E')dE'  dfi'  +  S  Ef(E) 

o 


Integration  over  energy  space  gives  angular  dose,  D  (r,  O),  due  to  the  angular 
neutron  energy  flux: 


°  ]  I  1 1  I  °  dE  dO 


+  (f,  n) 

where  Dj  is  the  dose  due  to  the  uncollided  angular  neutron  energy  flux  from  the 
source. 
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This  equation  can  be  written  as 


+  D^(r,  O) 


where  k(r',  O'  -•  r,  0)  is  the  dose  transport  kernel  and 


HP,  O'  -  *  O,  =  y  enE>K(E',  0',?^E,T,0). 


For  the  duct  problem,  our  concern  is  not  with  the  tronsport  medium  but  with 
reflection  from  the  walls  of  the  duct.  In  this  case,  the  dose  transport  equation  can, 
by  analogy,  be  written  as 


D  (r,  O) 


=  I  I 

O '  Area 


D(r',  0’)x(r',  O' 


r,  0)dAd0' 


+ 


(r,  O) 


(4) 


where  integration  is  carried  out  over  the  reflecting_^surface ,  A  (the  floor,  ceiling, 
and  walls  of  the  duct).  In  this  case,  the  kernel  x  (r'.  O'  "*  r,  O)  has  a  different 
meaning.  It  is  the  probability  that  the  angular  energy  flux  which  corresponds  to 
D  (r'.  O')  is  reflected  at  differential  area  dA  such  that  it  becomes  0(7,  O). 

Dj  (7  fl)  is  the  uncollided  dose  contribution  from  the  source. 

Consider  the  difference  between  Equation  3  and  Equation  4.  In  Equation  3 
the  quontity  of  interest  is  the  neutron  flux  as  a  function  of  angle,  position,  and 
velocity.  Therefore,  the  entire  range  of  neutron  energy  is  important.  However,  in 
Equation  4  the  important  quantity  is  the  neutron  dose  integrated  over  energy. 

For  the  same  value  of  n  (E,  r,  Q)  per  unit  energy  interval,  En  (E,  r,  fj)  always 
has  a  relatively  higher  value  in  the  high  energy  range.  Furthermore,  the  response 
function  for  neutron  dose,  f  (E),  has  higher  values  for  higher  neutron  energies. 
Therefore,  as  far  as  the  total  neutron  dose  is  concerned,  the  lower  energy  neutron 
contribution  is  much  less  than  the  higher  energy  neutron  contribution  if  there  is 
roughly  the  same  neutron  density  at  all  energies. 
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Now,  during  the  slowing-down  process  of  a  neutron,  while  it  is  streaming 
through  the  duct  by  reflection  from  the  walls,  the  more  reflections  there  are  from 
the  wall,  the  more  the  neutron  loses  its  energy,  until  it  is  slowed  down  to  thermal 
energy.  The  result  of  this  slowing-down  process  can  be  expressed  in  simple  mathe¬ 
matical  form. 

Let  1.  be  the  neutron  dose  contributed  by  the  ith  reflection: 

I  - 

lCT,  o)=  j*  J  I.  _^(;^,  n)x  (r",  O' -  T,  o)dAdn' 

Solid  Area 
angle 

i  (r,  O)  =  D  f,  O) 

o  s 

Since  neutrons  contributing  to  1;  have  been  scattered  at  least  i  times,  and 
since  each  scatter  results  in  a  degradation  in  energy, 

l.(T,  O)  >  O)  (5) 


Next,  consider  the  geometrical  factor  for  neutron  reflection  which  affects  the 
total  dose  In  a  duct.  For  simplicity,  consider  a  straight  cylindrical  duct  as  shown 
in  the  following  diagram. 


The  source  is  placed  on  the  centerline  at  the  entrance  to  the  duct.  The  detector  is 
at  r.  The  total  reflection  area  is  277’RL,  where  R  is  the  radius  of  the  duct,  and  L  is 
the  total  length  of  the  duct.  Then  the  dose  atTis 


(6) 


0(^,0)  =  y  J  D(?^,  n')H  o)dAdn'  +D^(^;n) 

Solid  Area 
angle 


Letting  t  be  the  axial  distance  from  source  to  detector,  Equation  6  can  be  written  as 
L  =  .C 

D(on)=  j  y  D(r‘^,  5')x  (i^,  O' ft)dn'27rRd^, 

Solid  L  =  0 
angle 

'  '  1.  - 

■  I  i'  «")h  (I^,  0"  -’o  0)d0".27rRd-t  +  0.  (^  0)  ::  (7) 

•  -.Solid 

;  ■  .■,..;angle: 

Now,  for  certain  values  of  -t,  the  . second  terrn  of, Equation  7  is  smaller  than, 
the  first  terms  for,  two  reasons.  (1)  The  distances  between  the  source  and  the  detec¬ 
tor  are'greater,  (2)  At  larger  distances,  higher  orders  of  reflection  ore  more  likely-  - 
relative  to  the  number  of  low-order  reflections. 

In  view  of  the  above  considerations,  the  first  approximation  to  the  solution 
of  Equation  4  Can  be  written  as 

ft)  ':w 


The  second  approximation  is 


^2^^  O)  J  J  ft’)x  (i^.  O'  -7,  0)dAdn'  +  ft)  (9) 

Solid  Area 
angle 
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The  third  approximation  is 


03(7,0)=  J  j  02(7,  o'))t(7,  O' -7  ft)dAdo' +  0^(7  n)  (lo) 

Solid  Area 
angle 

The  solution  of  Equation  4  can  then  be  written  as 


n  : 

0(7  0)  ^  0  (7  n)  =  y  ,1.(7  o)  (11) 

n  Z_i  '  ,  . 

i  =  o  -  . 


where  1^  (r,  O)  corresponds  to  0^  (r,  OV 

iiifa  <ai9lk  aJhk 

Since  Ij  (r,  O)  >  l-^  ^(r,  O),  as  previously  discussed,  the  series  approximt^ion 
of  the  solution  of  Equation  4  must  converge  under  oil  circumstances  because-  (r,  0) 
Is-bounded  by .ihe  source  strength  for  all  n.  . 

Equation  4  can  be  applied  directly  to  the  twaHegged  rectangular  concrete 
duct  problem.  The  dose  transpCrt  kernel  is 


a  (E^,  0  ,  0)  cos  0^ 
00  0 


where 


=  distance  from  source  to.  the  reflecting  area,  dA 
R2  =  distance  from  dA  to  the  detector 
cos  6^  =  cosine  of  the  polar  angle  from  the  source  dA 


a(E  ,0,0)  =  differential  dose  albedo 
o  o 
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In  the  second  leg  of  the  two-legged  duct,  the  source  function  of  Equation  4 
drops  out  if  the  detector  in  the  second  leg  cannot  see  the  source  directly.  Then 


D(^  ft)  =  y  l.(^  ft) 


■  ■  *•* 

Now,  for  the  singly  reflected  neutron  dose,  I.  (r,  ft)  carl  be  written  Cs 


j.  ■  /*  ^  -  D  cos  0 

I,  (r,  ft)  =  J  j  -^Y~r  ©o' 

.  -  -  Solid:  Area  ^ 
angle-  ' 


where  Dq  is  the  . dose  at  unit  distonce  from  the 
:  i :  “  Tht  doubly  ref Iscfed  neutron  dose  Is 


source  in  qir. 


.  (12) 


p  .  n  \/cOs.0 

j (r/  ft)-;=  r j-;:-:::  J  1^  (r'y  ft')  -  a  (E^/  0g>:  d)  dA  dft ' 

-"Solid -Area  .  .-  Hh  '  )  ■" 

— •"  . angle-  ■  ■  -  -  - • ..  .  ..-.r 


(13) 


where 


■=  distance  from  the  first  reflecting  area  to  the  second  reflecting 
area,  dA 


R2  =  distance  from  dA  to  the  detector 


cos6^  =  cosine  of  the  polar  angle  from  the  first  reflecting  area  to  dA 


a(E^,  6^,  0)  =  differentia!  dose  albedo  of  singly  reflected  neutrons 
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In  the  same  way,  the  neutron  dose  from  the  higher  orders  of  reflection  can  be 
obtained. 

The  general  formula,  Equation  4,  is  based  on  neutron  reflection.  However, 
there  are  a  certain  number  of  neutrons  which  penetrate  the  comer  lip  of  the  duct 
and  are  reflected  to  the  detector.  Also,  corner-lip  scattering  of  neutrons  contributes 
to  the  total  dose  in  the  second  leg  of  the  duct.  Therefore,  these  factors  have  to  be 
taken  into  account  for  the  t^a^dose  in  a  duct  with  a  bend. 

The  calculation  of  D  (r,  O)  was  carried  out  by  the  IBM-1620  computer.  The 
details  of  the  calculation  are  discussed  in  the  following  sections. 

1.  Calcuiation  of  Primary  Reflection 

Equation  12  is  equivalent  to  the  LeDoux-Chil ton  formula  for  calculating  the 
■  contribution  of  the  primary  reflection  to  the  gammd-ray  dose  in  a  two-legged  rec¬ 
tangular  duct. 7  In  their  calculation,  the  primary  reflecting  area  was  divided  as 
shown  in  Figure  4,  and  R],  R2,  Sq,  and  0  were  considered  constant.  For  this  calcu¬ 
lation,  the  primary  reflecting  areas  were  divided  into  five  sub-areas  as  shown  in 
Figure  5,  the  source  was  considered  to  be  at  the  mouth  of  the  duct,  and  the  detector 
was  token  to  be  placed  on  the  centerline  of  the  second  leg,.  The  detector  position 
coordinate  wos  (x^j,  y^,  o).  Each  sub-area  was  divided  into  m  increments.  The  dose 
rate  coritributed  by  the  reflection  from  each  sub-area,  AA,  was  calculated  using  the 
equation 


D.  (r,  ft) 


D  a.  cos  6  £iA, 
o  j  o,  J 
J  J 


J  J 


(14) 


The  dose  rate  due  to  primary  reflection  was  obtained  by 


D^(i‘ft)  =  y  0.(0  ft) 
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2.  Calculation  of  Secondary  Reflection 


Equation  i3  was  evaluated  in  the  same  way  as  Equation  12  was  evaluated  for 
primary  reflection.  The  secondary  reflecting  area  was  divided  as  shown  in  Figure  6. 
The  contribution  of  the  secondary  reflection  was  obtained  by  taking  combinations  of 
the  scattering  areas;  for  example,  the  combination  of  the  source  +  1  +  6  +  D,  or 
of  the  source  +  2  +  5  +  D.  This  time  the  scattering  area  was  not  divided  into 
Incremental  areas,  but  the  parameters  were  obtained  for  the  center  of  each  scattering 
area.  This  method  is  very  similar  to  that  of  Chopmon  for  calculating  gamma-ray 
dose  attenuation  in  two-legged  ducts. ^ 

The  working  equation  for  this  secondary  reflection  is 


D,  (7,  ft)  = 


and 


DjS  ft) 


3.  Corner  Inscattering  and  Transmi: 

The  corner  inscattering  and  transmission  contribution  from  the  corner  lip  of 
the  duct  is  discussed  in  great  detail  in  Reference  7  for  the  gamma-ray  cose.  For  the 
neutron  case,  some  modification  is  necessary  because  of  the  physical  nature  of  the 
problem.  The  corner  transmission  is  treated  the  same  way  as  it  was  for  the  gamma - 
ray  case.  The  treatment  of  corner  inscattering,  however,  is  based  on  the  assumption 
that  the  scottering  of  neutrons  in  concrete  is  isotropic  in  the  laboratory  system.  In 
this  calculation,  only  the  primary  effect  is  considered  for  corner  transmission;  and 
for  comer  inscattering  the  primary  and  secondary  effects  are  taken  into  consideration. 
For  the  secondary  corner  inscattering  effect,  only  four  special  cases  which  seem  to 
be  important  are  considered.  They  are  shown  in  Figure  7. 
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Figure  6.  Division  of  secondary  reflecfing  areas. 
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7.  Geometry  for  multiple  corner  Inscattering  with  sea; 
surface  in  the  first  leg  (Cases  1  and  2)  and  in  the  s( 
leg  (Cases  3  and  4).  Areas  2  and  5  In  Case  1  are  e 
ceil  ing  or  floor. 


CALCULATION  BY  MONTE  CARLO  METHOD,  USING  ADONIS  CODE® 

One  approach  to  the  solution  to  the  problem  of  neutron  strsaming  through  a 
duct  is  the  Monte  Carlo  method.  The  ADONIS  code  solves  the  transport  equation 
by  the  Monte  Carlo  method  for  neutrons  traveling  through  a  configuration  consisting 
of  a  finite  number  of  rectangular  parallelepiped  regions.  ADONIS  calculates  the 
flux,  F,  in  each  region  as  follows: 


F(r,j) 


S 


G 

where  S(r,  j)  =  (l/E  -  E  .)[W(r)/Gn]  V  S  (r,j) 

J  J  “  '  Z-j  9 

9=1 

track  length  of  the  gth  group  of  neutrons  in  region  r,  with  energy  E 
satisfying  E.  ,  ^  E  s;  E. 

■  J-l  J 

volume  of  region  r 
G  =  total  number  of  Monte  CorIb  groups 
n  =  number  of  neutron  histories  per  group 
W  (r)  =  weighting  factor  for  region  r 

The  flux,  as  defined  above,  is  the  total  track  length  per  unit  volume  per  unit 
energy  interval.  This  gives  the  flux  averaged  over  the  volume  of  the  region.  In 
order  to  get  the  flux  at  any  desired  point,  the  volume  of  the  region  must  be  suffi¬ 
ciently  small  so  that  the  flux  averaged,  over  the  volume  can  represent  the  flux  at 
the  midpoint  of  the  volume.  If  the  volume  of  the  region  is  small,  the  statistics  are 
usually  very  poor  unless  the  case  history  number  is  tremendously  high.  In  order  to 
overcome  this  difficulty,  it  is  assumed  that  the  flux  over  the  cross-sectional  area 
can  be  reasonably  well  represented  at  the  center  of  that  area. 

ADONIS  calculations  based  on  this  assumption  are  compared  with  albedo 
model  calculations  and  experimental  measurements  following  the  next  section. 


and  S  (r,  j)  = 
9 

V  = 
.  r 
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EXPERIMENTAL  MEASUREMENTS  OF  FAST  NEUTRON  DOSE  DISTRIBUTIONS  IN 
DUCTS 


Some  preliminary  experiments  have  been  carried  out  at  NCEL  by  Doty. 

Using  a  neutron  generator,  neutron  sources  with  energies  of  14  Mev  and  2.5  Mev 
were  obtained  by  T  (d,  n)  and  D(d,  n)  reactions.  The  dose  rates  were  measured  by  a 
£  12-inch  spherical  dosimeter  in  a  3  x  3-foot  concrete  duct.  All  the  values  measured 

^  by  the  dosimeter  were  normalized  by  a  monitor  value  obtained  with  a  paraffin- 

I  covered  BF3  counter  at  on  arbitrary  position.  The  results  so  obtained  are  reproduced 

I  in  Figures  8  and  9. 

I  In  this  experiment,  only  one  duct  size  was  used.  Therefore,  the  shape 

I  dependency  was  not  studied,  but  the  dose  distribution  as  a  function  of  the  distance 

5  from  the  intersection  of  the  two  legs  was  studied. 

I  The  dose  rate  fell  off  approximately  as  the  inverse  square  of  the  distance 

*  along  the  axis  of  the  firsr  leg  of  the  duct.  In  the  second  leg,  the  rate  fell  off  more 

I  rapidly  —  approximately  as  the  inverse  third  or  fourth  power  of  the  distance  along 

I  the  axis  of  the  leg.  Note  that  this  finding  Is  different  from  the  exponential  dose 

I  attenuation  that  would  be  expected  for  deep  penetration  through  a  homogeneous 

I  medium. 

I  Since  dose  rates  are  proportional  to  a  power  of  axial  distance,  it  seems  clear 

I  that  neutron  penetration  through  duct  walls  Is  small  compared  with  reflection  from 

j  the  walls.  The  above  argument  is  true  only  when  the  thickness  of  the  walls  is 

I  greater  than  a  mean  free  path  of  the  initial  neutrons. 

I  . 

r 

I  COMPARISON  BETWEEN  EXPERIMENT,  MONTE  CARLO  CALCULATIONS,  AND 

I  ALBEDO  CALCULATIONS 

i  Experimental  results  for  neutron  dose  were  obtained  from  the  spherical  dosimeter 

I  in  units  of  rem.  But  the  dose  calculated  by  the  albedo  model  is  given  in  rod  units. 

i  Therefore,  rad  dose  was  converted  to  rem  dose  to  allow  direct  comparison  of  theoret- 

I  ical  results  with  experimental  findings,  in  this  procedure,  the  average  energies  of 

I  the  singly  reflected  neutron  and  doubly  reflected  neutron  were  obtained  by  averaging 

J  Allen's  Monte  Carlo  data^  over  the  reflecting  angle  weighted  by  the  number  of 

i  neutrons  in  each  reflecting  angular  sector.  The  energies  thus  obtained  were  used  to 

^  convert  rad  dose  calculated  by  the  albedo  model  to  rem  dose  according  to  the  spher- 

t  ical  dosimeter  response  (which  is  close  to  being  tissue  equivalent). 9 

t  In  the  same  way,  the  results  of  ADONIS  Monte  Carlo  calculations  were 

I  converted  from  neutron  flux  to  dose  according  to  the  spherical  dosimeter  response 

!  function. 
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1  10  100 
Centerline  Distance  From  Source  (It) 

Figure  8.  Experimental  measurements  of  dose  for  T  (d,  n)  reaction 


Counts/Minute 


1  10  100 
Centerline  Distance  From  Source  (ft) 


Figure  9.  Experimenfal  measurements  of  dose  for  D  (d,  n)  reaction 
in  duct. 
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In  comparing  the  ADONIS  results  to  others,  in  order  to  evaluate  the 
representative  value  in  a  given  region,  the  following  technique  is  applied: 

Let  (^/^o)/i^QONIS  ~  ratio  of  dose  at  detector  to  dose  at  unit 

distance  from  the  source  in  air  obtained  by  ADONIS 
in  a  given  region 

(D/Dc)  =  results  of  the  same  ratio  calculated  by  the  albedo 

model  at  point  x,  which  is  in  a  region  corresponding 
to  the  region  specified  by  ADONIS 

The  mean  value  of  D/Dq  within  a  region  extending  from  xj  to  X|  +  ]  and  across 
the  entire  cross  section  of  the  duct  is,  by  the  Mean  Value  Theorem, 


I 


This  mean  value  for  D/Dq  can  be  obtained  by  graphical  integration  of  the 
values  for  (D/Dq)j^  obtained  from  albedo  colculations. 

Comparisons  con  now  be  made  between  D/Dq  and  (D/Dq),(^[5q|vj|5  for  the 
region  which  extends  from  Xj  to  x-^  The  comparison  is  made  at  the  jsoint  x" such 
that 


(f )  =  (f-) 

o  _  o 

X 


The  comparison  is  shown  in  Figure  10  for  an  1 1-inch  square  duct,  it  is  seen 
that  the  albedo  model  calculations  give  somewhat  lower  values  than  do  the 
Monte  Carlo  calculations.  This  result  is  expected  since  orders  of  reflection  higher 
than  the  second  are  neglected  in  the  albedo  model  calculations. 
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Figure  10.  Comparison  of  albedo  model  and  ADONIS  Monie  Carlo 
calculations  of  T  (d,  n)  reaction  in  a  1 1  x  1 1-inch 
two-legged  duct. 
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The  experimental  values*  and  the  albedo  model  calculations  for  the  T  (d,  n) 
reaction  (14-Mev  neutron  source)  are  shown  in  Figure  11.  The  geometry  of  the  duct 
used  for  the  calculations  and  the  experiment  was  not  exactly  the  same.  In  the  experi¬ 
ment,  the  source  was  placed  1  foot  outside  the  entrance  of  a  15-foot  first  leg;  in 
the  calculations,  the  source  was  placed  at  the  entrance  of  a  16-foot  first  leg.  In 
this  comparison,  it  is  assumed  that  the  geometrical  difference  of  1  foot  in  the  first 
leg  would  not  significantly  affect  the  calculated  results.  As  seen  in  this  comparison, 
the  calculated  values  are  slightly  lower  than  the  experimental  values.  Again,  this 
is  expected  since  the  higher  orders  of  reflection  are  neglected  in  the  albedo  model 
calculations. 

The  experimental  values  and  the  albedo  model  calculations  for  the  D  (d,  n) 
reaction  (2.5-Mev  neutron  source)  are  shown  in  Figure  12.  The  ADONIS  Monte  Carlo 
calculation  required  25  hours  of  IBM-7090  computer  time,  while  the  albedo  calcula¬ 
tion  was  performed  in  45  minutes  on  the  IBM-1620. 


FINDINGS  AND  CONCLUSIONS 

1.  Calculations  of  fast  ndutron  streaming  by  means  of  the  semiempirical  formula 
developed  in  this  study,  based  on  the  olbedo  concept,  are  in  close  agreement  with 
calculations  by  the  ADONIS  Monte  Carlo  technique. 

2.  Calculation  by  the  semiempirical  formula  took  about  one-thirtieth  the  computer 
time  required  for  the  Monte  Carlo  method. 

3.  Agreement  between  the  results  of  experiments  and  the  theoretical  calculations 

was  within  .15  percent.  --  _  . 

4.  For  the  purpose  of  obtaining  engineering  design  criteria,  it  is  sufficiently  accurate 
to  use  the  values  of  the  dose  distributio.ns.  in  the  second  leg  of  the  duct  which  are 
obtained  by  albedo  model  calculations  using  second-oider  approximations. 


FUTURE  PLANS 

This  was  a  preliminary  study  on  this  subject.  In  order  to  have  more  complete 
information,  it  is  necessary  to  study  the  problem  for  a  variety  of  duct  dimensions. 

In  the  near  future,  an  experimentol  study  of  corner  effects  and  the  effect  of  duct 
shape  will  be  compared  with  the  albedo  model  calculations.  Also  it  is  planned  to 
further  compare  the  three  independent  approaches  to  this  problem:  the  experimental 
studies,  albedo  model  calculations,  and  ADONIS  Monte  Carlo  calculations. 

In  order  to  obtain  better  results  with  the  albedo  model  calculations,  it  may  be 
necessary  to  calculate  the  higher  order  reflection  contributions  to  the  total  dose. 
However,  such  calculations  are  expected  to  be  very  difficult  and  to  require  lengthy 
computer  computations.  In  fact,  the  more  accurate  high-order  calculations  may  well 
be  prohibitively  expensive. 


D/D 


Figure  12,  Comparison  of  albedo  model  calculation  and  experimental 
measurements  of  D  (d,  n)  reaction  in  a  3  x  3-foot 
two-legged  duct  for  a  2.5-Mev  neutron  source. 
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Appendix 

COMPARISON  OF  RESULTS  OF  ONE-PARAMETER  FORMULA, 
SEMIEMPIRICAL  FORMULA,  AND  MONTE  CARLO  DATA 


Dr.  L.  B.  Gardner  and  Mr.  A,  J .  Mettler  of  NCEL 
furnished  the  results  of  the  Monte  Carlo  calculations. 
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13.  ABSTRACT  .  -■..r,.','- 

Asa  part  of  the  Laboratory's  fundamental  shielding  studies^  for  personnel  shelters/  fast 
neutron  dose  rates  are  colculated  in  the  second  leg  of  on  air  duct  through  concrete  for  neutron 
energies  of  14  Mev  and  2.5  Mev.  The  calculotional  technique  is  based  on  the  albedo  concept. 
Dose  rates  are  also  calculated  by  a  Monte  Corlo  technique,  and  the  results  obtained  by  the  two 
theoretical  methods  ore  compared  with  each  other  and  with  experimental  measurement's. 

Comparison  shows  very  good  agreement  among  these  three  independent  determinations. 
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